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We consider the Hubbard model on the infinite-dimensional Bethe lattice and construct a sys- 
tematic series of self-consistent approximations to the one-particle Green's function, G^ n \uj), n = 
2,3,... . The first n— 1 equations of motion are exactly fullfilled by G^ (uj) and the n'th equation of 
motion is decoupled following a simple set of decoupling rules. G^ 2 \lu) corresponds to the Hubbard- 
III approximation. We present analytic and numerical results for the Mott-Hubbard transition at 
half filling for n = 2, 3, 4. 
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C3 1 INTRODUCTION 

Si 

The Hubbard-modcl on the infinite-dimensional, half-filled Bethe lattice in the paramagnetic state has been consid- 
I ' ered in the last years more and more as the standard model for the Mott-Hubbard transition. One of the reasons for 
this interest is that the study of interacting Fermions in the limit of infinite spatial dimensions Q leads to considerable 
technical advantages, as the self-energy becomes stricly local. The Mott-Hubbard transition has been studied recently 
\ by Monte-Carlo ||, a self-consistent weak-coupling theory j^], rigorous j| and self-consistent exact diagonalization 

o 

ON ■ Alternatively, Hubbard H has considered an equation of motion approach to the Hubbard model. Unfortunatly his 
i approach can not be systematically improved in any finite dimensions, since the resulting self-consitency equations 
would be numerically intracable, due to the involved summations over momenta. Here we point out that, uniquely 
to the infinite-dimensional Bethe lattice, the equation of motion approach may be systematically carried on. We 
develop a simple decoupling scheme applicable to equations of motion in any orders and resulting in a self-consistency 
equation for the one-particle Green's function. We present results for the Green's function is second, third and fourth 



x 



studies Ira], with considerable different results, see 



order, with the second order corresponding to the Hubbard-III solution |fj. 

DEFINITIONS 

At half filling, n = 1, the chemical potential [i = U/2. For general fillings we set /i = U/2 + A/x and find for the 
grand-canonical Hubbard Hamiltonian, K = H — /j,N, 

K = t J2 W + tV2 1><-1) - A M J>, (1) 

<i 1 j> 1 a i i 

where the symbol < i,j > denotes pairs of nearest neighbour sites on a Bethe- lattice with coordination number z. 
The scaling t = t/^/z yields a non-trivial limiting behaviour in the limit z — > oo. In Eq. (Q) we have made use of 
cj „C; „, h, ; = h { ^ +n i I in Eq. (|l]) and of some of the following operator definitions: 
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where the q CT and the c ilT are creation/destruction operators for electrons with spin a =|,| on lattice sites i. The 
following operator identities hold: 



(W + ii,i, CT )/2 



For the equation of motion we will make use of following operator commutator relations: 

f r- 



tj,i,—ai j } — (7 



u ii3i,i,-o 



t i 



U i > 3j,i,-a C i,a 
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and similar one's. Given operators A and -B, the Matsubara Green's function ((A;B)) is defined as 



dr(-l) < T t 1(t)B(0) > e l "" T , 



with r being the imaginary time and the brackets < . . . > denoting the thermodynamic expectation value. The 
retarded Green's function is obtained via the analytic continuation iuj n — > uj+iS. Here we are interested in constructing 
systematic approximations to the onsite one-particle Green's function G(iuj n ) = ((c CT ; a )). In the following we will 
define sites 1, 1', 1", ... to be different n.n. sites of the central site, 0, and 2, 2', 2", ... to be different n.n.n. sites 
of the central site. 



EQUATIONS OF MOTION 

For any imaginary-time-dependent operator, A{t) in the Heisenberg picture, the equation of motion is given by 

jU(r) = [K,A], 

where, in the grand-canonical formulation, K = H — fiN is given by Eq. (Q). Note that we have set fj, = U/2 + A/z. 
Some examples of specific equation of motion in r-space are 
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where j is a n.n. site of i and where we have used definitions Eq. (|J) and the operator identity Eq. In Fourier-space 
the first equation of motion is 

{w n + A(i) «c 0iO .;cJ ;CT )) = l + 2t((c liff ;cJ iCT )) +C7/2((rf _ CT c 0iO .;4, CT )), 
where the site I is any of the z equivalent n.n. site of the central site 0. In the limit z — *■ oo the decoupling 

«c 1)CT ;qU> - ti^;^))^;^)) (13) 

is exact since all diagrams contributing to the propagation of the particle from site 1 to site 1 have only a vanishing 
probability 1/z — > to visit the central site. Using that ((c 1 a ; c\ a )) = ((c CT ; cj CT )) = G(ioj n ) in the paramagnetic 
state and that zt 2 — t 2 we rewrite the equation of motion of first and second order as 

[w n + Afi-PG(iu> n )]G(w n ) = l + [//2((d _ CT c 0i(T ;4 )(7 >) (14) 

(iw„ + A/x)«4- CT c 0jCT ;4 )CT )) = 2An ,_ o + U/2G(w n ) (15) 

where < c\ _ CT c _ a — c _ CT cJ CT > = 2rto,-cr — 1 = 2A?io,-o- and where we have measured no,-<r = 1/2 + Ano,- ff with 
respect to half-filling. The second-order equation of motion, Eq. ([l5]) , generates two new Green's function. In third 
order we have the equation of motion for _ a c a ; cj (T )), 

(iw n + An) «j'i,o,- CT c 0i0 .;4 l0 -» = ((c{_ a c _ a - 4-<A,- CT » + (O'i.o,-^^,^; 4 )ff » 

- Z H(k,0.-a £ ,a^lj) + Zt((i lv C 0a ,cl )) 

(16) 

+ 2*((n 1) _ (7 c 0i<7 ;4 )ff )) - 2t ((n ,-a £ o.a^oJ) 

where we have used various commutator- relations, in particular Eq. (||). We can simplify Eq. @ in infinite-dimensions 
where the decoupling 

«%,-<A<7;4,<t» < > ((c , CT ;4 )O -)) ( 17 ) 

becomes exact. Again < % _ ff > = 1/2 + An^-^ and using Eq. (Q), 2n _ CT — 1 = c? we can rewrite 

2*<K- CT e 0jCT ;4,<r>>- 2t«VA«r5 0) = 2 *An 1 ,_ CT ((c 0i<T ;4 iCT » - * (18) 
in Eq. (fi6|). In third order we have with the equation of motion for ({d _ a c 1 a ; c\ a )) a second equation of motion, 

(iuj n + Ap) {{d Q c x -,cl )) = zt{{d _ a c 2 ]cl j) -2zt((j vo c l a ;cl a )) 

(19) 

+ i«VA*;4,„» + ^/2<<vA-A^ a U>- 



A total of six new Greens's functions are generated in third order by Eq. ( [l6| ) and Eq. (|H|), ({d _ a d x _ cr c 1 a ; 4 a }} 
and «rf _ CT c 2)CT ;cJ )CT )), by Eq. @ alone, )) by Eq. 

( |l6| ) alone, and (Oio-oA' o-'^o.o-)) ^y both Eq. ( |l6| ) and Eq. (|l9|). We present the equation of motion for these six 
Green's functions in Appendix A. Now we specialize on half-filling, where Afi — 0. For an overview we rewrite here 
the first two equation of motion, Eq. ( |l4| ) and Eq. ([l5]) : 



[ic; n - PG(iw„)] G(iu n ) = l + U/2((d _ a c l7 ;cl j) (20) 

iuj n ((d Q _ a c a ;cl a )) = U/2G(w n ) - 2zt ((j 1Q _ ff c 0i(T ; + ^ «4,-<A >0 -;4,<7»- ( 21 ) 

Noting that there is no current flowing in equilibrium, < cj _ CT c _ a — cj _ cr c 1 _ CT > = and using Eq. ( [i"8| ) with 
Arii.-cr = the first third-order equation of motion, Eq. (|l6|), takes the form 

iUn(Gi,o ,-*Co.*;4>.*)) = *t{{j xfi c v \c\ )) - ((4.o -^0,^4 )) + zt ((i l v c ;4 )0 -» , x 

(22) 

For completeness, we rewrite also the second third-order equation of motion, Eq. (|l6|), for the half-filled case, 

((4o,-<A >er ;4 = rf ((4-a a 2.,;cL)) ~ 2 2* {(j ll0 )) 

(23) 

+ *<<4,-A<^ 6 U> + ^/2«rfo,- CT rfi,- CT c 1>CT ;4 iCT )). 

Particle-hole symmetry at half- filling implies that the real part of G(u>) is an odd function of frequency. Consequently 
Eq. (|2(]) is a even function of frequency, Eq. ( pl| ) an odd function of frequency and Eq. ( |22| ) and Eq. (53) even functions 
of frequencies. The absence of a (real) constant term on the right-hand side of Eq. (|Tj) is therefor forced by symmetry 
but the absence of constant terms on the right-hand side of Eq. (E2L) and Eq. (J23|) is accidental. 



APPROXIMATIONS 

In order to obtain a solution for G{iLo n ) we will decouple the equations of motion. We present here a simple, 
systematic decoupling scheme, which may be applied at any order of the equation of motion. We generalize the 
decoupling for ((c lff ; 4 a)) presented in Eq. (|l3] ) by 

(Qi,o,-Joy^o,a)) -> -tG(iu n ) ((d c ;cl )) 

' ' „' (24) 



which may be motivated by considering Eq. ([22|) and Eq. ([23|). In contrast to Eq. (|13j), which is exact in the limit 
of infinite dimensions, Eq. ( p4| ) is an approximation only, valid in the paramagnetic state. In a state with staggered 
antiferromagnetism the replacemment G(iLu n ) — > G- a (iuj n ) should be made on the right-hand side of Eq. ([24]). Using 
Eq. ( pil| ) and Eq. (^0|) we then obtain 

[iu n - ?G(iu n )] G(iu n ) = 1 + - G(iw n ). (25) 



Comparing Eq. (25) with the Dyson equation of the infinite-dimensional Bethe- lattice, which has the form 
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[iw n + fj, — - PG(iw„)] G(iw n ) = 1, (26) 

we find 

£< a >(iw„)-M = r- (27) 

iu> n — 3 t 2 G(iu! n ) 

Eq. (|27]) is the well-know result of the Hubbard-III approximation ||. Here we use the index ^ to indicate that for 
this approximation to the self-energy the first equation of motion has been retained exactly and that the equation 
of motion of second order has been decoupled. In order to obtain £( 3 )(iu; n ) we have to consider the equation of 
motion of third order, Eq. J22] ) and Eq. (|23|). For the decoupling of the terms ~ t on the respective right-hand sides 
straightforward generalization of Eq. (J2J) can be used. For the terms ~ U/2 we proposed the dcoupling scheme 



(28) 



valid in the half- filled case, where we have been motivated by Eq. flAip and Eq. Eq. ( p2| ) and Eq. ( Bjj ) then 

become 



+ 3?G(iw n )«ii 

+ 3?G(t« n ) (vV<A,^4, CT » + ( s (^n) - aO (vVA^L))- 
We replace £(«<;„) — > £( 3 )(iw„) in Eq. ( |29| ) and obtain for Eq. ( pl| ) 

3p 

iw„((do c ;cj » = {7/2G(iw n ) + - nr 2r , r — ; T^IWT- — ^ V i(d c ; cj )). 



(29) 



Finally, using Eq. (E6J) and Eq. (20) we obtain 



S (3) (-n) - /I = , (30) 

ia; n - 3 t 2 G(iu n ) - (£( 3 ) (io; n ) - /i) 
with /z = {7/2 at half-filling. In Appendix A we derive the expression for £( 4 )(iw n ), see Eq. (All). Considering the 



self-energy expanded in powers of (U/2) 2 and of t 2 , we see that T,^(iu! n ) is exact up to t 2 and all powers of (U/2) 2 . 
Similary Y,^(iu> n ) is exact up to t 4 and T,^ (iuj n ) up to t 6 . 



RESULTS 



The solution for the retarded Green's function, G(uj), is given by Eq. (E6J) together with the respective self-energy 



see Eq. (|27|), Eq. ( p0[ ) and Eq. (All). For small {7's the solution is metallic, for large {7's insulating. The transition 
point can be calculated from a small-u; expansion, by noting that in the insulating state the Laurent-series of the 

self-energy starts by a 1/w divergence [Q: S(w) ~ a/u + From Dysons" s equation, Eq. (p6|), it follows than that 

that G(w) ~ —iv /a + . . .. We find then that 
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= = \{U 2 -\2t 2 



a 



(4) _ j U 2 - 160 t^/U 2 (31) 



I" 



where has been calculated in Appendix A (see Eq. ( A12 )). It can be shown [Q, considering the spectral repre- 
sentation of G(oj), that the coefficient a > 0. The U c for the Mott-Hubbard transition is then given by the a = 
condition in Eq. (|l]). We then find [7 C (2) = C/ c (3) = VT2t and J7 C (4) = \/l60 1. From now on we use the scaling 
t = 1/a/2. We have 

UP = C/ c (3) = a/6 ~ 2.4495 



[/ c (4) = ^40 - 2.5149. 



(32) 



The ?7 C given in Eq. ( |32j ) has been denoted j8| t/ c i, being the critical interation strength at with the insulating solution 
becomes stable. A U C 2 has also been defined || as the critical interaction strength at with the metallic solution becomes 
unstable, i.e. when ImG(0) = 0. We will discuss further below that the equation of motion solutions always yield 
£/M = t/ c ( " } ee U { c n) as given by Eq. ©, for n = 2, 3, 4. 

In Fig. |l](a) we have plotted a/ 3 ^ and a' 4 ' as given by Eq. ([H]). For comparision we have included in Fig. [l] the 
results obtained by an exact diagonalization study Q, where the symbol S(l) denotes the so-called Hubbard star ^ 
and the symbol S(2) the so-called star of the stars. In general, the symbol S(n) denotes clusters which are truncated 
Bethe-lattices of order n. For instance, S(l) is the cluster containing a central site with its z n.n. sites. As the 
S(n) are finite clusters, no true Mott-Hubbard transition is observed, only a crossover. Nevertheless, good agreement 
between the S(2) cluster and the equation of motion results if found for U > 3, indicating that U c < 3. 

In Fig. |l|(b) we have plotted the Z-factor, as defined by the seond term in the Laurent expansion for the self-energy, 
S(o;) ~ a/uj + (1 — 1 / Z) lj + . . ., which we calculate in Appendix B. In a Fermi-liquid state (a = 0) the Z-factor would 
have the meaning of the inverse effective mass, Z — m/m*, but in the insulating state (a > 0) considered in Fig. [j](b) 
the Z-factor is just a parameter. Z vanishes at U c for all three equation of motion solutions in the same fashion. To 
see this we note that «W ~ U- ui n) (for U ~ ut ] < ui n) and n = 2, 3, 4, see Eq. ©) and Z^ - [a^] 2 for small 
a M (see Eq. @ and Eq. (IbJ)). We then obtain 
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for n = 2, 3, 4 and C/ C (n) given by Eq. ©. 

In Fig. U and Fig. || we present the equation of motion results for the density of states, —Im (uj + iS), for a small 
8 = 0.0001. The results for U = 1 and U = 2 are in the metallic state and for U = 3 and U = 4 in the insulating state. 
The magnitude of the gap changes little with n = 1, 2, 3, but the shapes of the Hubbard bands change somewhat and 
side bands appear at higher energies. 

In Fig. ^ we present the results for the density of states at the Fermi level, -JmG'"'(0), and the self-energy at 
the Fermi level, -ImS'"'(0), in the metallic state, as a function of U. Both quantities are related via Y,( n \Q) = 
fx - G<")(0)/2 - 1/G(™)(0), compare Eq. @, with i?eS(™)(0) = ^. Analytically we find 



-ImG(2)(0) = 

V. tZl (33) 



! (^ + 3/2(172)") • 



G 



In a Fermi-liquid state JmE(O) = and — JmG(O) = \/2 [(To), independent of interaction strength, U. As we see form 
Fig. U -/mGW(O) < \/2 and the —Im £(™)(0) > 0. The equation of motion solutions do not describe a Fermi liquid. 
Indeed, in the equation of motion solutions the divergence 



lim -JmS<"'(0) ; 1 -» oo (34) 

drives the Mott-Hubbard transition. Here we have derived Eq. from Eq. ( p3| ) for n — 2, 3 and verified it for n = 4 
from a small G(0)-expansion of Eq. ( All ). We note that the metallic solution becomes unstable, as seen from Eq. 
([m|), at exactly the same Uc at which the insulating solution becomes stable, compare Eq. (|32]). In an alternative 
approach to the Mott-Hubbard transition on the Bethe lattice in infinite dimensions || it has been proposed that the 
metallic solution might remain stable up to a much higher U C 2 ~ 4.2 — 4.7. 



CONCLUSIONS 



We have considered the equation of motion approach to the Hubbard model on the infinite-dimensional Bethe-lattice 
and shown that it is possible to evaluate higher-order Green's function by a simple decoupling scheme. We presented 
analytic and numerical results for the half-filled case in the paramagnetice sector. We found the critical U c ~ 2.5 
for the Mott-Hubbard transition to change only little with the order of the approximation. The metallic state is a 
non-Fermi-liquid, with a finite imaginary part of the self-energy at the Fermi-level, diverging at the Mott-Hubbard 
transition. The metallic state does not shown, surprisingly, any tendency to become more Fermi-liquid like with 
increasing order of approximation. 

The author would like to thank G. Hiilsenbeck for discussions and a carefull reading of the manuscript. This work 
was supported by the Minister fur Wissenschaft und Forschung des Landes Nordrhein-Westfalen. 



APPENDIX A: 



Here we present the equations of motion of fourth order for the six new Green's functions created by the third-order 



equations of motion, Eq. (16) and Eq. (|19|). The first one is 

(iw n + A/ 1 )((4,- ff rf v / v ;cJ iff )) = U/2((d Q ^c lr7 ;clj) + 2tAn l ,- <r {{d 0i _ ir c 0itr ;^ itr )) 

+ zt ((vA-A*; 4,a>) - 2 zt (O'i'A- 4,<r» (A1) 

where we have used the exact decoupling 

{(d _A-<A<r; 4,<7» -> < d l.-<x > «4) -*\<t\ 4,a)) 

with < d 1 _ a > — 2Ani._ cr . Note, that the n.n. site 2 occuring in the last term of the right-hand side of Eq. ( |Al| ) in 
j 2 i _ a could in finite dimensions also be the central site, 0, but not in infinite dimensions. The second equation of 
motion in fourth order is 
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(iu n + A/i) {(d h J li0 _ a c y, c f 0a )) = < d h J lfi _„ > - U/2 {Oi,o,-cr%^ d oJ) 

+ ^(KAa-A',,;^))- 2 ^{Ck,i,Ji,o-<T% a -A,^)) ^ 

where we have used, besides others, the commutator relation [t 1 _ a ,t 1 _ a ] = 0. The third equation of motion in 
fourth order is 



(A3) 



(iw n + A/i) ((d ,- CT c 2)CT ; cJ, CT )) = t {{d _ a c h(T ; cj, ff )) + zt ({d^_ a c 3 al cJ >CT )) 

- 2 «5'i',o,-<A^ cL» + U/2 A- A»5 clj). 
The fourth equation of motion in fourth order is 

(iu n + A/i) (0 V>O _ a c ljCr ; cl j) = t {Q Vfi _ a c^ a ; 4 >CT » + z* ((?i>, ,-<A CT ; 4,<r» 

+ 2tAn 1 /,_ ff «c ii(T ;4 1< r»-*«4),-A«r;4 >ff )> 

- ((*2', 0> -^ 1)CT ; 4,J> + ** (&M",-A«r; 4 )CT » ^ A4 ) 

- (7/2 «4 A',o,-A<x; 4,«r» + ^2 «4A ,0,-A*; a U>' 
+ ^/2<(7V,o,-A-A«r;4, ff »> 

where we used the exact relation 

similar to one valid in Eq. ( |l6| ) and Eq. (|l^) . The fifth equation of motion in fourth order is 

(^ n + A/i)((t 11 , _ (T C 0i(T ;C^ CT )> = <i 1 ,l',- CT >+*(0i,O,- CT C Ol(7 ;4, ( r))- *((jo,l',-a 6 0, CT ;4,<r)> 

(A5) 

- 17/2 «d lt4r t l , 1 ,,_ or a , <T ; 4,«r» + U/2{{d v J lv _ a c^-cl a )), 

which might be further simplified due to t _ ff c 0j(T ; cj )0 .» - * <<5'o,i',-<A -; 4,<r>) = 2 * ((JiA-^o,^ 4,<r» and 
-^((^.l'.o-Co.o-^o.J) + «*((3'i,2'- -Co, -;4, -)) = - 2 ^((j2,i', ff CQ,<7;4 )0 -))- Finally, the sixth equation of motion in 
fourth order is 

(iuj n + A/i) «£,,0,-<A<7; 4,<r)) = < *2,0,- CT > - 1 (Q\,0,-<r%,<,) 4,<r)) + z * ((?2,1' -<A^ 4,<r)> 

- ^((J 3 ,0,<tC 0!(t ;4,<t» + zt ((h,0,-^V,a'A,o)) ^ 

- U/2{{d 2 J 2fi .c .-,clj). 



Due to particle-hole symmetry the real part of Eq. (A2) is an odd functions of frequency at half filling, where A/i = 0. 



Therefore the constant term of the right-hand side of Eq. (IA2), < d 1 rT i 1 _ a >, has to vanish at half-filling, and it 



does (compare the definition of d lcr , Eq. (||)). Note, that the difference v _ <T c Q CT ; cj — ((i 2 _ (T c CT ; cj g.)) is 
created in the third-order equation of motion, Eq. (p2p, and that therefore the difference of the constant terms on the 



right-hand side of Eq. (A5) and Eq. (A6) respectively, namely < t 1 > — < t 2 _ CT > has to vanish at half-filling, 
and it does. 

Now we derive Y,^ (iu} n ) at half-filling. Generalizing the decoupling schemes Eq. ( pi[ ) and Eq. (|2^) we define 

U/2 

(A7) 

b(iu n ) 



a(iu) n ) 



iuj n — 5 1 G(iiv n ) 
t 



and find 



iuj n — (£(iu>) — fi) — 3 P G(iu> n ) 

a(iw n ) ({d ;cj _)) 



((4 ,-cA,-«A,«r; 4,<r» = 

((4', CT *l',0,-<r 6 )( r; 4,<r» = ((il'.O-^O.tr! 4,<r)) 

((ii'A-o-ci^.'Co.a)) = & ( ia v> ( ((i'i',0 .-o-co.o-; 4,<7» - ((4),-<A, CT ; 4,<r») 

- CT C , CT ;4 l( r)) - ((Vo,-<A^4,<t>) = 26 ( iw n)((5'l'0,- CT e 0, O -;4, O -))- 



Here we have introduced another decoupling rule, namely the self-energy H(iu) n ) comes in definitions Eq. (A7) only 



(A9) 



with the prefactor zero or minus one. We substitute Eq. (A8) into Eq. ( |22| ) and Eq. (|23j). We find 

2B(iu n ) {Q v cl tCT )) + A{iu n ) «d ,-<A, CT ;4,<T» = 4 «4),-<r e o 1<7 ; 4 )0 .», 

with B(ioj n ) = ztb(iu> n ) and A(wj n ) = icj n — U/2a(iuj n ) — 3ztb(icJ n ), with a(iuj n ) and b(ioj n ) given by Eq. (|A7| ) 
Inverting Eq. (A9) we obtain 

((il',0,-cr^0,<ri ^0,<t)) = 

((<4-<A, ff ; 4,<7» = 



A(^„) + B(iuj n ) 



t A 2 (iuj n )-2B 2 (iuj n ) 
A(iuj n ) + 2B(im n ) 



(0o, 



-CT U 0,ff' 



4,.)) 



(A10) 



'A 2 (^„)-2^(^„) 
Finnally, using Eq. (pH), Eq. (|2^) and Eq. ( p6| ) we obtain 

(C//2) 2 



((^O.-aCQ^^L)) 



E (4) (iw„) - /i 



• , _ r 2 2 A(w») +35(M n ) ' 
A\iu n )-2B*(iu n ) 



(All) 



The critical interaction strength, Uc 4 \ at zero temperature may by be obtained from Eq. ( All ) by considering the 
insulating state. For small real frequencies cu the the self energy starts like E(u>) = a/uj + O(oj) and the Green's 

function like G{w) = -uj/a + 0(uj 3 ), which leads to B(w) ujP/a, A(lu) (a/oj)(U/2) 2 /(a + 5t 2 ) and £( 4 )(w) ~ 

{U/2) 2 /(u n -2t 2 /A{uj)). One then obtains 

U 2 /A-A{)t A /U 2 



a 



(4) 



l + 8t 2 /U 2 



(A12) 



Since a > Eq. ( |Al^ ) is valid only for U > Ur> = </l60i. For the usual scaling t = l/y/2 the critical U { c > = </40 
2.5149, which compares to uP = C/ c (3) = \/6 ~ 2.4495. (compare Eq. @). 



9 



APPENDIX B: 



Here we consider the insulating state, in which the Laurent expansion of the retarded self-energy starts like 

EH = g + /i + (l-^)w + ... 



where we have used Eq. (Eq) for self-consistency. In order to calculate Z we need 



(Bl) 



1 



u + /i-S(w)-3f G(w) 
1 



w-5P G(w) 



= w/a - (1 - l/Z)u 3 /a 2 + ... 

= -u/a - {1/Z + 3t 2 /a)oj 3 /a 2 + 

= " _ 5 P 1 /^ + *!/? • ■ 

(a + 5?)w (a + 5P) 2 



Using Eq. (E7h, Eq. (BO) together with Eq. (66) and Eq. (B2) one finds 



(B2) 



For the evaluation of Z^ we need 



2A{u) + 3_B(^) 



Finnally obtain 



Q 



(4) 



[a( 2 >] 2 /[([//2)V 2 ' +3i 4 ] 
[a( 3 )] 2 /[(C//2)V 3 ' +9t 4 ] 



a + 5P 



a(f//2) 
(a + 5P) 2 



p 



a 2 (C//2) 



7 (1/Z + P/a)^ 3 



,(4) 



a 2 (J7/2) 4 

(aW+5?) 2 



(2 + 9P/a)w 3 + 



[l + 8t 2 /U 2 ] = (U/2) 2 + P 



(U/2) 4 



2 + 9P/a( 4 ) 



10 P 



([//2) 2 a< 4 » 



from which we may obtain Z^ 4 \ whith Eq. ( A12) for 



(B3) 



(B4) 
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FIG. 1. (a) The coeficient a and (b) the coefficient Z of the self-energy, E(o>) ~ a/uj + U/2 + (1 — 1/Z)u> + . . ., in the 
insulating state at half-filling, as a function of U. Shown are the results for E' 2 ^(cj) (dotted line), £^(u;) (dashed line) and 
E (4) (u>) (solid line). For comparision, the corresponding results obtaind by exactly diagonalizing the Hubbard star, S(l) and 
the star of the stars, S(2), are given. 



FIG. 2. The density of states, -ImG {2) {u)) (dotted line), -JmG (3) (u) (dashed line) and -lmG (4) (w) (solid line) in the 
metallic state at half-filling for (a) U = 1 and (b) U = 2. 



FIG. 3. The density of states, ~ImG {2) {uj) (dotted line), -JmG (3) (ui) (dashed line) and -JraG (4) (u) (solid line) in the 
insulating state at half-filling for (a) U = 3 and (b) U = 4. 



FIG. 4. (a) The density of states at the Fermi level in the half-filled metallic state, -7mG (2) (0) (dotted line), -7mG (3) (0) 
(dashed line) and -JmG (4) (0) (solid line) as a function of U and (b) the self-energy at the Fermi level, -7mE' 2 '(0) (dotted 
line), -7mE (3) (0) (dashed line) and -JmE (4) (0) (solid line). 
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